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Abstract:  Formal  progressing  wave  expansions  are  applied  to  problems 

Involving  one-dlmenslonal  wave  propagation  through  lnhomogenous  elastic 
media.  Expansions  for  the  stress  and  particle  velocity  are  obtained  In 
addition  to  the  expansion  for  the  particle  displacement  which  Is  a special 
case  of  previous  results.  Although  one-dlmenslonal  problems  could  be 
solved  with  the  previously  reported  asymptotic  methods,  it  Is  more  con- 

Iff. 

venient  to  use  the  expansion  In  terms  of  the  stresses  to  evaluate  the 
expansion  coefficients.  The  procedure  Is  illustrated  by  solving  several 
problems  In  layered  and  nonlayered  Inhomogeneous  media  where  the  compresslonal 
wave  speed  Is  subject  to  power  and  exponential  variations  with  distance. 
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1.  IimODOCTIOM 


In  recent  yaara , tha  tachniquaa  of  geometrical  optics  have  boon 

gonoraliced  to  treat  tho  electromagnetic  wave  equations,  the  acoustic 

wavs  aquations,  tha  elastic  wave  equations,  and  tha  viscoelastic  wave 

equations,  chiefly  by  workers  at  the  Courant  Institute  of  Kathaoatlcal 

Sciences.  The  method  was  originally  applied  by  Lunabura  [13]  for  tine 
periodic  waves  in  the  form  of  an  asyaptotlc  expansion  in  inverse  powers 

of  the  frequency  to  show  that  the  leading  tarn  of  such  an  expansion  for 

tha  electromagnetic  wave  equations  is,  in  fact,  the  geometrical  optica 

solution.  It  was  shown  by  several  authors  [3,  7,  8,  10,  11]  that 

subsequent  terms  give  corrections  to  the  geometrical  theory.  The  method 

was  generalised  by  Frladlandar  [3]  to  treat  progressing  wave  forms. 

In  this  form  the  method  is  a convenient  tool  in  studying  the  propagation 

of  discontinuities. 

Karal  and  Keller  [5]  extended  the  method  to  treat  general  progressing 

waves  in  Inhomogeneous  elastic  media.  In  their  analysis,  the  formulation 

was  in  tarns  of  displacements  and  displacement  potentials.  For  certain 

problems,  it  would  be  more  convenient  to  have  the  formulation  in  terns 

of  velocities  or  stresses.  Bence,  for  the  one-dimensional  case,  these 
results  are  presented  here,  and  the  method  is  demonstrated  by  making  use 
of  the  stress  formulation. 
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It  should  be  pointed  out  that,  while  the  procedure  is  formal,  the 
results  for  the  leading  term  are  consistent -with  those  obtained  by  the 
theory  of  weak  solutions  and  other  more  fundamental  means.  For  a more 
complete  discussion  of  these  methods  see  Keller  (6]  and  references 
contained  therein. 

After  the  technique  is  developed,  one-dimensional  waves  that  propagate 

into  an  inhomogeneous  elastic  medium  of  constant  density  but  varying 

compresslonal  wave  speed  are  studied.  It  is  expected  that  these  results 
may  be  useful  for  studying  a surface  generated  wave  <fts  it  .propagates 

downward  into  the  near  surface  regions  of  the  earth's  crust.  Belated 

problems,  solved  by  different  techniques,  have  been  recently  treated  by 

Llndholm  and  Doehl  [12],  Whittier  [IS],  Payton  [13].  Other  related 

references  are  contained  in  these  papers. 

The  method  for  treating  interfaces  between  different  inhomogeneous 
media  is  then  presented.  As  an  example  case,  the  method  is  applied  to 
treat  an  inhomogeneous  layer  overlying  a homogeneous  half  space,  where 
a step  pressure  is  applied  to  the  free  surfece. 


2.  FORMULATION 

The  one-dimensional  equation  of  motion  for  a linear  elastic 
continuum,  assuming  small  strains,  is 


o^  - P(»)ott  • P<*)*t 


where  o,  p , u,  and  v = are  stress  (positive  in  tension)^  density 
particle  displacement,  and  particle  velocity;  where  x and  t denote 


spatial  position  and  time;  and  where  a subscripted  variable  denotes 
partial  differentiation  with  respect  to  that  variable.  Hooke's  law 


for  elastic  media  is  given  by 


where  £ is  the  modulus  of  elasticity.  The  variables  p(x)  and  E(x) 


are  assumed  to  be  continuous,  differentiable  functions  of  x.  Taking 


the  partial  derivative  of  (2)  with  respect  to  t,  «e  obtain 


Hence,  the  equations  for  one-dimensional  elastic  wave  propagation  in 


terms  of  displacements,  velocities,  and  stresses  are  given  by  combining 


We  nov  assume  that  the  displacement,  velocity,  and  stress  solutions 


where  the  prime  denotes  total  differentiation  with  respect  to 
the  argument , In  this  manner,  all  of'  the  f 's  may  be  related  to 


f.  (called  the  waveform)  by  successive  integrations.  For  example,  if 


fn  is  the  Heaviside  step  function,  H(t),  then 


r 


i 


Note  that  fQ  vanishes  for  negative  argument,  i.e.,  in  front  of  the 

vavefront  whose  equation  is  given  by  t ■ S(x),  where  S(x)  is  called 

the  phase  function.  For  this  case,  note  that  the  coefficients  U V or 

n,  n 

PQ  are  the  jump  conditions  for  the  displacement,  velocity,  or  stress 
and  their  derivatives  across  the  wavefront.  It  is  assumed  that  U 

o, 

V , and  P are  identically  zero  for  n < 0. 
n n 

From  the  relationships  in  (2),  (3),  (7),  and  (9)  with  the 
definition  v • ut,  if  may  be  shown  that  the  coefficients  in  the  displace- 
ment, velocity,  and  stress  series  are  related  by 


Vi  • 

(12) 

e(K  - s'Vi>  • 

(13) 

(14) 

P 


3.  SOLUTION 


Consider  the  displacement  solution  first.  If  (7)  Is  substituted 
Into  (4),  then 


m 

l 


U(s’)2e  - p]  unfn_2 


[2ES'U  * + (ES"  + E'S')U  ] f , 
n n n-l 


+ [EU  ••  + e'U  '1  f } - 0 . 
n n n 


(15) 


Since  Un  ■ 0 for  n < 0,  the  summations  may  be  changed  so  that  each 

term  Is  In  the  form  of  a coefficient  times  f . Then,  because  f Is 

n n 

related  to  f^  as  Illustrated  by  (11)  and  because  of  the  arbitrariness 
of  fQ,  esch  coefficient  may  be  equated  to  zero;  hence 

[(S')2e  - p]  0 ..  - 2ES’U  ' - (ES"  + E'S')  U 
htjl  n n 

+ EU'*  + E'U’  . - 0 . (16) 

n— l n— l 

Setting  n " -1,  noting  that  U >0  for  n < 0,  and  requiring  without 

A 

loss  of  generality  that  UQ  ^ 0,  we  obtain 

(S')2  - p/E  - c’2  (17) 


7 


which  nay  be  recognized  as  the  eiconal  equation  of  gemoetrlcal  optics 
where  c(x)  Is  the  wave  speed.  Hence, 


S(x)  «St 


dx 


where  S = S(x)  Is  the  Initial  value  of  the  phase  function,  l.e.,  Its 
value  at  x ■ i.  In  (18),  the  plus  sign  Is  taken  for  waves  which  propagate 
In  the  positive  x direction  and  the  minus  sign  for  waves  which  propagate 
In  the  negative  x direction.  Substituting  (18)  into  (16),  we  obtain  a 
linear  first  order  ordinary  differential  equation  whose  solution  is 


U (x)  - U 
n n 


fp(x)E(x)  f , 

P (x)  c (x)  + J n-1 


p(x)c(x) 


where 


n-1  2pc 


£ <EU’n-l>  ‘ 


In  a similar  manner,  the  solution  for  the  coefficients  of  (8) 


and  (9)  are 


ejilslO.  df 

p(x)c(x) 


r — — 

f 

4.  EXAMPLE  PROBLEMS  FOR  INHOMOGENEOUS  HALF-SPACES 
Consider  the  half-space  x > 1,  defined  by  the  properties  E - E(x) 
and  p ■ p ■ constant,  initially  at  rest.  For  t > 0,  a pressure  is  applied 


to  the  surface  x ■ 1 ao  that  o(l,t)  is  known.  The  solution  will  first 
be  obtained  for  a step  pressure  input  and  Duhamel's  theorem  will  be 
applied  to  obtain  the  solution  for  a general  input  preasure  function. 


The  phase  function  S Is  then  determined  from  (18)  and  the  general  f 

n 

la  obtained  from  (11) . Since  the  density  will  be  constant  In  the 


example  cases,  from  (24)  reduces  to 


C 


n-1 


1 

2 


(27) 


Case  I.  c(x) 


- 2a 

c x ; p - p 


In  this  case,  E 


- -2  4a  . „ 
pcx  and  C 


n-1 


may  be  shown  that 


cx2a  P"..  From  (18),  It 
~~2  n-i 


S(x)  - x~2a  * 1-1 
c(l  - 2a) 


. 


- iog  X , 

c 


a - r 


(28) 


From  (23),  (25),  and  (27),  It  may  be  shown  by  induction  that  if  a i*  ■=■ 


a i y . 


(29) 


j-0 

Substituting  (29)  Into  (23)  and  simplifying,  we  obtain 


11 


..  


I 


Ajn  " 2 Aj-1,  n-1 


La  + 2a-l)  (1-1)  H-o  + 1 (2q-l)  1 l<j<n 


J (2a-l) 


.-£  f A [a  + (2e-l)(J-l)][-a  + 1(2a-l)1 

2 jil  AJ’l.n-I  j(2a-l) 


1 - 0,  n > 0, 


-1  , 

0 , 


J - n - 0 , 
j < 0 or  J > n , 


whore  ..  5 V If  “ * 2"*  it:  “I  again  ba  shown  by  induction  that 


K -V*  l Ain(l0«  *>J  . 
j-o  J 


(31) 


where 


AJn  “ 2 I(J  + X)  AJ  + 1,  n - 1 " 4j  AJ-1,  n-11  , 1 - i * n , 
• 0 , J i 0,  n > 0 or  J > n , 


-1,  J - n - 0 . 


(32) 


All  poaaibla  values  of  a have  bean  included  in  the  above  results. 
In  particular,  the  case  a • 1 is  of  Interest  because  it  aay  be  shown  f 


12 


J 


(30)  that  A__  - -1  and  A.  >0  fpr  all  other  values  of  j and  n.  Hence, 
uu  jn 

the  series  terminates  after  a single  term, 

o - -x  H [t  - (x-1)  /cx]  ,o-l,  (33) 

which  may  be  shown  to  be  the  exact  solution  for  this  case.  Note,  in 

this  case,  that  the  stress  is  unbounded  for  increasing  x,  and  that,  for 

ct  > 1,  a - -x  since  (x-l)/x  < 1 for  all  positive  x.  That  the  stress 

should  grow  without  bound  is  non-physical.  However,  we  note  that  the 
_ 2 

assumption  c - cx  is  also  non-physical  for  unbounded  x.  Hence,  if  the 

_ 2 

model  c * cx  is  realistic  for  x < L < • , it  is  expected  that  the 
solution  (30)  is  applicable  for  x S L < - . 

From  the  above  results,  it  is  seen  that  the  stress  amplitude  may 
either  increase  or  decrease  with  increasing  x depending  on  the  value  of 
a.  In  the  case  a - 0,  the  above  results  reduce  to  the  homogeneous 
elastic  case  where  the  stress  level  is  independent  of  the  distance 
traveled.  For  o > 0,  the  stress  level  Increases  and  for  a < 0,  the  stress 
level  decreases  with  increasing  x. 

Case  II.  c - c exp{2B(x-l)  } , p-p 

In  this  case,  E - p c^  exp  (4B(x-l)}  and  ■ j exp  {20(x-l)> 

P".  . From  (18),  it  may  be  shown  that 
n— i 
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(34) 


r 


S(x)  - 1 -Jgp 
2 Sc 


From  (23),  (25),  and  (27),  It  may  be  shown  by  induction  that 


PB  - exp  (B(x-l)}  l A exp  (26J(x-l)}  . (35) 

“ j-0  Jn 


Substituting  (35)  into  (23)  and  simplifying,  we  obtain 


, . m-i>* ; t * 

J»  *J  J-l,  n-l  • 


1 < i < *»  , 

• ^ I-*’*- 


-1  » 
0 


j - n - 0 , 

J < 0 or  J > n . (36) 


Again,  the  above  solution  reduces  to  the  limiting  homogeneous  elastic 
caaa  for  6-0.  As  x increases,  the  stress  grows  if  6 > 0 and  attenuates 
if  S < 0. 

In  the  above  examples,  the  response  was  determined  for  a step 
input  compressive  stress  at  x - 1.  The  stress  solution  has  the  form 


14 


H [t  - S<»>] 


wh«r«  Pn(x)  Is  given  by  (29),  (31),  or  (35)  and  S(x)  Is  givsn  by  (28) 
or  (34).  In  general,  tbs  solution  for  an  arbitrary  boundary  strsss 
o(l, t)  nay  bs  determined  from  Duhamal's  theorem 


ft  3oH(t) 

J 


o(l,  t - t)  dt. 


5.  REFLECTION  AND  TRANSMISSION  OF  WAVES  AT  AN  INTERFACE 
Consider  a stress  wave,  traveling  In  the  direction  of  Increasing 
x In  a medium  "a,"  Incident  at  an  Interface,  x - x,  between  medium 
"a"  and  medium  "b."  Let  the  compress lonel  wave  speed  and  density  of 
the  two  media  be  denoted  by  c#,  c^  o#,  and  pfc.  It  Is  assumed  that 
a wave  Is  transmitted  and  a wave  Is  reflected.  Due  to  the  linearity 
of  the  equations,  the  stress  field  can  be  obtained  by  superposition. 
The  boundary  conditions  at  the  interface  require  that  o - o.  and 

A D 

».  ■ t.  . Assume  that 

m D 


”*  ‘Jo  '«•  '*  (t  • Sl»  \l  '«  <*  - V • 


°b  "Jo  *»  " ‘ *»>  ' 


Where  the  subscripts,  1*1,  2,  3,  denote  the  incident,  reflected,  and 

transmitted  waves  respectively.  Since  the  form  for  all  of  the  waves  is 

that  represented  in  the  previous  sections,  the  solution  is  known  providing 

the  initial  values,  S , P , V.  , are  known.  It  is  assumed  that  the 

i in  m 

incoming  wave  has  been  defined  by  the  methods  of  the  previous  section. 
Hence,  S,  and  P^n  are  known. 

From  (14)  and  (17)  it  can  be  deduced  that 

pcvjn  - (-1)^  (PJn  - ,n-!>  1 ■ • <«> 


where  F=  pc2  v’  . (42) 

J .n-1  j .n-1 

and  s’  • c-1  for  the  outgoing  waves  and  s'  ■ -c”1  for  Incoming  wave. 

Note  that  F.  . is  a function  of  the  coefficients  P.  where  k s n-1. 
j »n-l  ^ 

In  particular,  (41)  yields  the  veil  known  result  from  the  conservation  of 
momentum  across  a steadily  moving  discontinuity,  P^  - -pc | Q | . 


When  equations  (39)- (42)  are  combined  with  the  boundary  conditions. 


o ■ a.  and  v <■  v.  , we  find  that 

a d a d 


S1  " S2  " S3  * 


and 


P2n  “ P3n  ’ _Pln  » 


P2n  + rP3n  " Pln  * Fl,n-1  + P2,n-1  + rF3,n-l, 


where 


' " V*b  S 


Hence : 


(43) 


(44) 


(45) 


P2n  * Pln  - Fl,„-1  + P2,n-1  * * F3,n-1>  (1  + • 


P3n  " [2Pln  " F1»n_1  + F2»n_1  + r F3»n“l1  (1  + r>1 


(46) 


Thus,  since  and  P^q  are  known,  Che  initial  values  at  the  Interface 
(x  * x)  are  determined  and  the  procedures  of  the  preceding  section  may 
be  applied  to  obtain  the  formal  solution.  One  should  be  careful  to 
choose  the  proper  sign  in  (18)  , (20) , (22)  , and  (24)  so  that  the  direction 
of  propagation  is  accounted  for.  It  will  be  noted  that  the  above  results 
reduce  to  the  limiting  homogeneous  elastic  results  in  the  proper  way. 

One  interesting  result  from  (46)  is  that.  In  the  limiting  case 
of  r ■ 1,  there  is  no  reflected  wave  only  if  F,  , ■ F_  , ■ F,  , 

i im  |Q*X  J 

for  all  n > 1.  Thus,  even  though  the  impedance  is  continuous  at  x it  is 
still  possible  that  a wave  will  be  reflected  if  the  impedance  of  medium  b 
varies  in  a different  manner  than  that  of  medium  a.  Note  that,  in  the 
case  cfl  ■ c^  , the  reflected  wave  has  no  discontinuity  since 
P20  ” °*  Father,  its  amplitude  increases  from  zero  at  the  front  in  a 
manner  which  depends  on  the  functions  representing  the  impedance  in 
media  a and  b.  This  general  result  was  pointed  out  for  a particular  case 
by  Payton  [13]. 

- - 2 

As  an  example  case,  consider  that  for  l<x<x^,  p>p  and  c * cx 

_ 2 

where  c • c at  x • 1.  For  x > x^  let  c " cx^  and  p - p . Thus,  the 
Interface  is  between  an  inhomogeneous  layer  overlying  a homogeneous 
elastic  half  space.  Further,  the  Impedance  is  continuous  at  the 

* - • 
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Interface.  Consider  the  result  when  s step  pressure  Is  applied  at  the 
surface  x - 1.  The  Incident  stress  wave  Is  given  by  (33).  The  stress 
field  is  given  by  the  first  of  (39)  In  the  layer  and  by  the  second  of 
(39)  In  the  half  space.  This  solution  will  apply  until  the  time  that 
the  reflected  wave  reaches  the  surface  x - 1.  At  x - the  phase 
functions  are 


xl~1 

S1  " S2  " S3  * - 
C*1 


and  from  (18) 


x.  (1  + x)  -2x  xf  -2x.  + x 

s2  " : >3 --4-4 

c Xj  X e x^ 


It  may  be  shown  from  (41)  and  (46)  that,  since  P1Q  - -x  and  Pj.  - 0 
for  n > 0,  - 0 for  n > 0;  f »*  0 for  n > 1;  P2n  - ?3n  - ^n-l^2 


for  n > 1;  and 


P9_  ■ -(-cx./2)x  , nil 


, n £ 0 


p3n  ■ -(-cx1/2)nx1  , n > 1 

■ -x,  , n ■ 0 


, n < 0 


Thus,  the  solution  Is  determined  by  substitution  of  (48)  snd  (49) 
Into  (39).  As  was  pointed  out  previously  this  solution  Is  applicable 
only  until  the  time  that  the  reflected  wave  strikes  the  surface  x • 1. 

At  this  point  the  procedure  may  be  repeated  for  that  boundary  and  so  on. 


6.  SUMMARY 

A formal  procedure  for  obtaining  solutions  for  one-dlmenslonal 
wave  propagation  In  Inhomogeneous  elastic  media  has  been  presented 
and  Applied  to  two  general  cases.  There  is  no  proof  that  the 
resulting  expansions  are  convergent  to  the  exact  solution;  however, 
in  all  cases  where  exact  solutions  are  known  and  where  the  technique 
has  been  applied,  it  has  been  found  that  the  correct  asymptotic 
solutions  have  been  given.  In  one  of  the  example  cases  treated  here, 
the  series  terminated  and  a closed-form  exact  solution  was  obtained. 

The  method  of  treating  a single  interface  between  different 
Inhomogeneous  media  was  presented  in  detail.  Multiple  Interfaces  may 
be  treated  In  exactly  the  same  way.  One  interesting  result  was  that 
reflacted  waves  may  be  generated  in  the  case  of  en  incident  wave  at 
the  point  where  two  different  Inhomogeneous  media  are  joined  even 
though  the  Impedances  at  the  Interface  are  matched. 


The  method,  while  tedious,  is  direct  and  Involves  only  ordinary 
differentiation.  Integration,  and  algebra.  The  tens  of  the  series 
are  evaluated  recursively  once  the  Initial  values  of  the  coefficients 


\ 


1 


and  the  phase  function  (the  barred  variables)  are  known.  These  Initial 
values  are  detenlned  froa  application  of  the  boundary  conditions 
in  a straight  forward  Banner. 
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